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Abstract. The Euler equations governing a relativistic perfect fluid are put into 
symmetric hyperbolic form with dependent variables the fluid's specific entropy 
plus a generalized velocity vector equal to the fluid's unit relativistic velocity 
vector divided pointwise by the fluid's specific enthalpy. 

Introduction. Frauendiener pQ and Walton [2j independently derived covariant 
symmetric hyperbolic Euler equations for isentropic relativistic perfect fluids. The 
purpose of this paper is to derive covariant symmetric hyperbolic Euler equations for 
the more general, nonisentropic relativistic perfect fluids. 

The dependent variables of a relativistic perfect fluid include a unit velocity 
vector field U a plus two thermodynamic scalar fields which I will choose to be the 
fluid's specific enthalpy h and specific entropy s. Following Frauendiener and Walton, 
U a and h are merged into a generalized velocity vector field J a = U a /h. Then, 
the Euler equations for the new dependent variables J a and s are derived from the 
assumptions that energy-momentum and material are locally conserved by the fluid's 
motions. These Euler equations are not symmetric hyperbolic, but they imply the 
fluid's specific entropy s is necessarily a first integral of J a . The latter is used to rewrite 
the Euler equations in an equivalent symmetric hyperbolic form. Having symmetric 
hyperbolic Euler equations for a relativistic perfect fluid is important for two reasons. 
First, existence and uniqueness of solutions to the Euler equations are immediately 
guaranteed by symmetric hyperbolicity. Indeed, proving existence and uniqueness 
for solutions of the usual non-symmetric hyperbolic Euler equations is difficult; see 
Lichnerowicz [3] for examples. And, second, to numerically solve the Euler equations, 
symmetric hyperbolic equations are preferred. 

In this paper Greek letters are used to denote tensor indices and units are used 
in which the speed of light = I . 

Thermodynamics. In this paper always assume a relativistic perfect fluid obeys 
an equation of state p = p(h,s), defining the fluid's pressure p to be a function of 
the fluid's specific enthalpy h and specific entropy s. The fundamental relation of 
thermodynamics, 

dp = n dh — nT ds, (1) 

then informs us that the fluid's material density n and absolute temperature T can 
be derived from the equation of state by the equations: 

-i < 2 » 
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and 

T = - l - 9 £- (3) 
n us 

And, from the definition of the fluid's energy density e = nh — p, the scalar field v 
derived by the equation: 

-2 / ue \ de/dh hdn h d 2 p 
\dp J dp/dh n dh ndh 2 ' 
is the speed of sound in the fluid. 

Constitutive Functions. Let g a p be the spacetime's metric tensor and assume it 
has signature — h ++. And, following Walton |2] and Frauendiener P, let J a be a 
vector field tangent to the timclikc fluid flow and normalized to define the value of the 
fluid's specific enthalpy by the equation 

h=(-g af3 J a J (3 r 1/2 . (5) 

Then the relativistic velocity vector field U a , energy- momentum tensor field T£, and 
material current vector field N a of a relativistic perfect fluid are the following functions 
of the fields J a and s: 

U a = hJ a , (6) 

T^ = nh 3 J a J^ +pgi, (7) 

N a = nhJ a . (8) 

And the tensor field 

q a /3 = g a p + h 2 J a Jp (9) 

has the useful algebraic property of projecting vector fields onto their component 
orthogonal to J": q^jP = 0. 

Some Calculus. Let V M denote the covariant derivative operator with symmetric 
connection defined by the Christoffel symbols of the metric tensor g a p, so V p,g a p = 0. 
Using equations (QJ, J5J), (@J, and (JSJ we can derive: 

V fl h = h 3 JpV fl J P , (iO) 

VflU = dh V>lh + 9I VaiS = nh2v ~ 2j P v ^ jP + ^" v m s > ( n ) 

V M j3 = nV^/i - nTV^s = nh 3 g^Jp\/ v J - nTV M s. (12) 

And, then from the definitions of the energy-momentum tensor field and the material 
current vector field we can compute the covariant divergences: 

V/^V.^JJ" + pg£) 

= (y fl n)h 3 J a J tl + 3nh 2 (V^h)J a J^ 1 

+ nh 3 {g a pj» + J a g^)V^ + ng£V - nTg^V^s 

= nh 3 {v' 2 h 2 J a J fj J^ + Z^JaJpJ* 1 + g afj J^ 

( dn \ 
+ J a 4 + &Jp)V vt J f> + i-^h 3 J Q J^ - nTg£) V p s 

= nh 3 (v- 2 h 2 J a J fJ J^ + q aP J» + J a q» + q£ J /3 )V Al J 13 

— h 3 J a J»-nTgZj\7^s, (13) 



and: 

V^iV = V p (n/iJ") 

= 7ih(v~ 2 h 2 JpJ" + h 2 JpJ" + g^V^J + ^hJ^V^s 

= nh(v- 2 h 2 J p J^ + gg)V„ J" + ^ WV„s. (14) 

Euler Equations. The Euler equations are derived from the assumption that the 
fluid's energy-momentum tensor field and material current vector field both have zero 
covariant divergence, so energy-momentum and material are each locally conserved. 
From equations l|13|) and Q14JI. the Euler equations are explicitly: 

nh 3 (v- 2 h 2 J a J J» + q aP J» + J a q% + ggJ^VpJ" 

dU -h 3 J a J» - nTgA V pS = 0, (15) 
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and: 



nh{v~ 2 h 2 J p J^ + q^)V^J P + ^/iJ^s = 0. (16) 

Equations (|1 5|) and i|16[l are not symmetric, and therefore not symmetric hyperbolic, 
because the coefficient of V M ,s in (|15|l and the coefficient of V M J@ in (|16(l are different. 

Symmetric Euler Equations. From the Euler equations (|15|) and it follows 
that the specific entropy s is a first integral of the vector field J a : 

J Q V M T£ + V^nhJ^) = -nTJ^V^s = 0, (17) 

so a relativistic perfect fluid flows adiabatically. Multiples of J^V^s can therefore be 
added to/subtracted from the Euler equations (|15[) and (|16fl and the new equations 
will be equivalent to the Euler equations, if the new equations also imply s is a first 
integral of J a .% In particular, subtracting \{dn/ds)h + (1 + v^ 2 )nT]h 2 J a J fJ, S7 
and [{dn/ds)h + (1 + v~ 2 )nT] J M V M s from equations l(T5|) and ifTCl) . respectively, 
and multiplying the resulting second equation by —T/h will symmetrize the Euler 
equations. Equations (| 1 51) become: 

nh 3 (v~ 2 h 2 J a J p J» + q aP J» + J a q% + q^V^J 13 

~nT(v- 2 h 2 J a J^ + <?£)V M s = 0; (18) 

and equation (|16f) becomes: 

- nT(v- 2 h 2 Jf 3 J» + gpV^J^ + nT 2 h- x (\ + v~ 2 ) J^V^s = 0. (19) 
Because the linear combination of equations (jlSJ) and (|19(l : 

J a ■ (Ella - hT~ l ■ (QU = -nTJ^V^s = 0, 

implies the specific entropy is a first integral of the vector field J a , the symmetric 
equations i|18|) and i|19fl are equivalent to the Euler equations (|15|l and (|16|) . Equations 
l|18fl and (|19f) will be referred to as the symmetric Euler equations. 

Isentropic Euler equations. A relativistic perfect fluid is isentropic if its specific 
entropy s = constant. According to equation (|17|l . in an isentropic perfect fluid local 

X Equivalent means fields J a and s which solve the new equations will also solve the Euler equations 
1151 and 1161 . and conversely. 
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material conservation is a consequence of local energy-momentum conservation 
11511 . An isentropic relativistic perfect is therefore governed by the reduced system of 
equations: 

nh 3 (v- 2 h 2 J a J J» + q ap J» + J a q£ + <£J /3 )V„J /3 = 0, (20) 
with only the field J a as a dependent variable. Frauendiener pp and Walton [2] 
independently derived equations l|2UI) and proved they are symmetric hyperbolic. 

Characteristic Analysis. Introduce an orthonormal vector field basis consisting of 
the unit fluid velocity vector field U a = hJ a plus three unit spacelike vector fields 
X a , Y a , and Z a . And let ^ = luU^ + K^U^ = 0, be an arbitrary vector field. 
There is no loss of generality by assuming the vector field = kX^, k > 0. Then, 
with respect to this basis and to the vector field £ M , the characteristic matrix for the 
symmetric Euler equations (|I8|I and l|19|) is the symmetric 5x5 matrix: 
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T 2 /i- 4 (l+v- 2 )w 
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(21) 



And the characteristic polynomial equation (dispersion relation) for equations (|f 8|l 
and (Jini is: 

det (A £ ) = -n 5 h 6 T 2 Lu 3 (v- 2 Lu 2 - n 2 ) = 0, (22) 

where det(— ) denotes the determinant operator for square matrices. Obviously, 
for each wave vector = kX^ the characteristic polynomial equation has only 
real roots in the frequency variable lo, if the speed of sound v is real valued: 
ll>i(k) = lo2{k) — u>3(k) — 0, ll>4(k) — vk — v^/kFk^, and ^(k) = —vk — —v^/kFk^. 
Therefore, according to the theory of partial differential equations (Courant & Hilbert 
0]): equations and are hyperbolic with wave modes propagating in the 

bicharacteristic directions: 

xJ = U a + g^^M =U * + ^Mr, 7=1,2, 3, 4, 5, (23) 
ok' 1 k 

tangent to three dimensional characteristic hypersurfaces with normal vector fields: 

# = wi(k)U>* + kX", 1=1,2, 3, 4, 5. (24) 
Such wave modes propagate at speeds Vi(k) = \u>j(k)\/k, 1 — 1,2, 3, 4, 5, with respect 
to observers comoving with the fluid at the relativistic velocity U a . The modes 
I = 1,2,3 are stationary with respect to the fluid; they represent flow along the 
fluid's streamlines. The modes 1 — 4,5 propagate at the speed v with respect to 
the fluid; they represent sound waves, because no other wave modes are possible in a 
perfect fluid, and justify interpreting v as the speed of sound in the fluid (Taub 0, 
Lichnerowicz 

Symmetric Hyperbolicity. As already noted, the characteristic matrix (|21ll of 
equations l|18|) and l|19|) is symmetric. Let W be a real valued, 5x1 matrix with 
transpose W T = (Wi, W 2 , W 3 , W 4 , W 5 ). Then equations (JTSJ and HHJ) will be 
symmetric hyperbolic with respect to a vector field £ M if the quadratic form defined 
by the characteristic matrix (|21|l : 

W T • A e • W = -nh 2 {[v- 2 (W! + Th- 2 W 5 ) 2 + W% + W% + Wj 

+ T 2 h- 4 Wi]uj 

+ 2{W X + Th- 2 W 5 )W 2 K}, (25) 
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is definite for all W 7^ (Courant & Hilbert 0], Friedrichs 

Suppose the vector field £ M = loU^ + kX^ satisfies the inequality [v^q^ — 
h^JuJuj^C = v 2 k 2 — uj 2 < 0. The latter is equivalent to > vk > 0, implying 
uj 7^ 0. Definiteness of the quadratic form (|25|) then follows because: 

-W T • A 6 ■ W/(ujnh 2 ) = v- 2 (W! + Th- 2 W 5 ) 2 + W| + W% + Wj 

+ T 2 h- 4 Wg + 2(Wi + Th- 2 W 5 )W 2 n/uj 

= (|Wi + Th~ 2 W 5 \/v- \W 2 \) 2 

+ 2\Wi + Th- 2 W 5 \\W 2 \/v + Wl + Wl 

+ T 2 hT 4 Wl + 2(Wi + Th~ 2 W 5 )W 2 K/uj, 

has a lower bound which is a sum of nonnegative summands with at least one summand 
always nonzero if W 7^ 0: 

-W T • A ? • W/(ujnh 2 ) > + Th~ 2 W 5 \/v - \W 2 \) 2 + W 3 2 + Wj 

+ T 2 h~ 4 Wl 

+ 2\Wi + Th- 2 W 5 \\W 2 \{l/v - k/\lj\). 

Here | • | denotes the absolute value of the enclosed quantity. 

Therefore: if the speed of sound v in a relativistic perfect fluid is real valued, then 
the symmetric Euler equations llfif) and li.9|) are symmetric hyperbolic with respect to 
any vector field £ M which satisfies the inequality: (w 2 ^^ — ft ^^J^f < 0. 

Interestingly, symmetric hyperbolicity of the Euler equations puts no limit on the 
speed of sound in a relativistic perfect fluid other than it be finite. But physics does: 
in any physically reasonable relativistic perfect fluid, the speed of sound should not 
exceed the speed of light: v < 1. 

The Cauchy Problem. For any scalar field r such that V^r 7^ 0, the equation 
r = constant defines a local foliation of spacetime by three dimensional hypersurfaces 
with normal vector field g^ u V v r. The Cauchy Problem for the symmetric Euler 
equations (|18|l and i|19fl asks: does there exist a unique set of fields J a and s which 
both satisfy equations l|18|) and l|19|) and take freely specified values J a = J Q | T= o 
and s = s| T= o on the initial value hypersurface t — of the foliation? From the 
works of Friedrichs 0, Courant & Hilbert 0], Fischer & Marsden [HUH], and others: if 
equations llty) and are symmetric hyperbolic with respect to the normal vector field 
£P = g^ v V ' v t of the foliation, and if the initial values J a and s have square integrable 
derivatives of order N > 4 on the initial value hypersurface, then the answer to the 
Cauchy Problem for equations ffflj) and $19\) is yes. 

A Final Note The symmetric Euler equations (|18|l and (|19(l for the fields J a and s 
are not uniquely symmetric hyperbolic. For any positive scalar field 9, one can add 
(9 - l)nT 2 h- 1 J fl V fJ ,s to the left hand side of J33J and obtain: 

- nT{ir 2 h 2 J p J^ + q^VftJ 13 + nT 2 hr 1 (9 + iT 2 )J M V^s = 0. (26) 

The system of equations ((18(1 and (|26[) is symmetric, equivalent to the Euler equations 
11511 and and symmetric hyperbolic with respect to any vector field £ M which 

satisfies the inequality: (v 2 q^ — h 2 J^Jv)^^ < 0. Fields J a and s which satisfy 
equations (|18J) and ((26(1 therefore also satisfy the equations ()15(l and 1(16(1 , or equations 
118( 1 and ((19(1 . and conversely, and so must be independent of choice for the arbitrary 
positive scalar field 9. Thus, there is no loss of generality by setting 9 = 1 and 
considering equations ((18(1 and l(19|) to be the symmetric hyperbolic Euler equations 
for relativistic perfect fluids. 
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